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Abstract. In this work we show that the systems of balance equations (balance systems) of 
continuum thermodynamics occupy a natural place in the variational bicomplcx formalism. 
We apply the vertical homotopy decomposition to get a local splitting (in a convenient do- 
main) of a general balance system as the sum of a Lagrangian part and a complemental "pure 
non-Lagrangian" balance system. In the case when derivatives of the dynamical fields do not 
enter the constitutive relations of the balance system, the "pure non-Lagrangian" systems 
coincide with the systems introduced by S. Godunov [Soviet Math. Dokl. 2 (1961), 947-948] 
and, later, asserted as the canonical hyperbolic form of balance systems in [Miillcr I., Rug- 
geri T., Rational extended thermodynamics, 2nd ed., Springer Tracts in Natural Philosophy, 
Vol. 37, Springer- Verlag, New York, 1998]. 

Key words: variational bicomplex; balance equations 
2010 Mathematics Subject Classification: 49Q99; 35Q80 



1 Introduction 

Systems of balance equations (balance systems) are a cornerstone of the mathematical models 
of continuum thermodynamics. That is why the usage of variational methods for the study of 
such systems was attempted in the corresponding domains of classical physics for a long time. 
For a review of the present state of this research field, we refer to the monographs |19l [20] and 
the references therein. These works use the conventional setting of variational calculus based on 
the action principle. 

Infinitesimal alternatives to the action principle were suggested by A. Poincare and E. Cartan 
in mechanics, by C. Caratheodory, H. Weyl and Th. Lepage in field theory and gradually evolved 
into a mathematical theory of a rare beauty, see, for instance, [H HI O [12] . Thus, it seems natural 
to investigate and to develop methods of infinitesimal variational calculus that can be used in 
the study of general balance systems. In the papers |13[ [T4"] it was proved that any balance 
system allows a natural realization by an appropriately modified Cartan form (see Section [3]) . 
Recently, I've observed that the balance systems present the largest class of systems of partial 
differential equations that can be realized using the forms of Lepage type on the infinite jet 
bundle J°°(7r) of the configurational bundle ir [15] , 

One of the most successful approaches in variational calculus is the method of variational 
bicomplex, see [TJ [T2] . In this paper we show that general systems of balance equations occupy 
a natural place in the framework of variational bicomplex. 

In Section [2] we introduce the necessary notions and notation. In Section [3] the balance 
systems of order k are defined. We recall the construction of the modified Poincare-Cartan 



*This paper is a contribution to the Special Issue "Symmetry, Separation, Super-integrability and Special 
Functions (S 4 )". The full collection is available at [http://www.emis.de/journals/SIGMA/S4.htmll 



2 



S. Preston 



forms for these systems, leading to the intrinsic formulation of the system of balance equations. 
An (n, l)-form Kq carrying all the constituents of a balance system is introduced. 

In Section|3]we show that any balance system is realized by the functional form corresponding 
to the (n, l)-form Kc- Higher versions of balance systems are introduced here as well. 

In Section [5] we apply the vertical homotopy decomposition to get a local splitting of a ge- 
neral (n, l)-form Kc into the Lagrangian part (corresponding to a function we call "quasi- 
Lagrangian") and a complemental "pure non-Lagrangian" balance system (K-splitting). Some 
properties of this decomposition are presented, including induced splitting of corresponding 
functional forms (F-splitting) . Finally, we obtain the condition on the form Kc to have zero 
Lagrangian part. 

In Section [6] we determine the conditions under which the Lagrangian part of the functional 
form for a balance systems of order zero (when the derivatives of dynamical fields do not enter 
the constitutive relation) vanishes. It is proved that in such a case, a "pure non-Lagrangian" 
system is of the type introduced by S. Godunov in 1961 [5] and later asserted by T. Rug- 
geri and his coauthors, as the canonical (symmetrical hyperbolic) form of balance systems in 
rational extended thermodynamics [11] . We calculate the "quasi-Lagrangian" for the balance 
systems of hyperelasticity and the 2-dimensional ideal plasticity and show that its input into 
the corresponding balance system reduces to the source terms in the balance laws. 

2 Settings and notations 

Throughout this paper ir : Y — > X is a (configurational) fibre bundle with a n-dimensional 
connected paracompact smooth C^-manifold X as the base and the total space Y, dim(Y) = 
n + m. The fiber of the bundle ir is a m-dimensional connected smooth manifold U. 

The base manifold X is endowed with a pseudo-Riemannian metric G. The volume form of 
metric G will be denoted by rj. Our considerations are mostly local and, as a result, we shall 
not deal with the boundary of X. 

We will be using fibred charts (W, x^,y l ) in a domain W C Y. Here (ir(W),x^) is the chart 
in X and y l are coordinates along the fibers. A local frame corresponding to the chart (W, x^,y l ) 
will be denoted by (d^ = <9 x m, <9j = d y i) (the shorter notation will be used in more cumbersome 
calculations) and the corresponding coframe (dx^,dy l ). 

We introduce the contracted forms rj^ = ig p ry. Below we will be using following relations for 
the forms r\ v without references (here and below \q = ln(-y/|G|)) 

dx v At]^ = 5^7], dr]^ = X G ^V- 

Section s : V — > Y of the bundle ir, defined in a domain V C X, corresponds to the collection 
{y l (x)} of classical fields. Usually these fields are components of tensor fields or tensor densities 
fields (sections of "natural bundles" [1]). 

For a manifold M we will denote by tm '■ T(M) — > M the tangent bundle of a manifold M 
and by r| 7 : T*(M) — > M its cotangent bundle. For a bundle Y ->■ X, denote by ir v : V(ir) ->• Y 
the subbundle of tangent bundle T(Y) consisting of the vertical vectors i.e. vectors £ € T y {Y) 
such that ir^y^ = 0. 

Denote by Vl r (M) — > M the bundle of exterior r-forms on a manifold M and by Q*(M) = 
®5^Lgf2 r (M) the differential algebra of exterior forms on the manifold M. 

Given a fiber bundle ir : Y — > X denote by J k (ir) the k-jet bundle of sections of the bun- 
dle 7r [HI [16]. For a section s : U — >• Y, denote by j k s : U — > J k {ir) the A;- jet of section s. 

Denote by irtr '■ J k ( 7T ) ~ > ^i^), k > r > the natural projections between the jet bundles 
of different order and by ir^ : J k (ir) —> X the projection to the base manifold X. 
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Projection mappings TTk(k-i) '■ «^ fc ( 7r ) J k 1 ( 7r ) form the tower of fc-jet bundles 

> J k (ir) -> J fc "V) -> (2.1) 

A bundle 7Tfc(fc_i) : J k (ft) — > J fc ~ 1 (vr) has the structure of an afftne bundles modeled on the 
vector bundle S k T*(X) ®^ Jfe -i (7r) F(tt) -»• J fc_1 (7r). 

Denote by J°°(7r) the infinite jet bundle of bundle ir - the inverse limit of the projective 
sequence (|2.ip [91 \12\ I16j . Space J°°(7r) is endowed with the structure of inverse limit of 
differentiable manifolds with the natural sheaves of vector fields, differential forms etc., making 
the projections ir^k : J°°{^) —> J (tt) the smooth surjections. See [5J [21 [16] for more about 
structure and properties of A;-jet bundles. 

For a mutliindex / = {i\,...,i n }, i& G N denote by 9 7 the differential operator d 1 f = 
cf 1 ! • • • 9*" in C°°(X). To every fibred chart (W, x^,y l ) in Y there corresponds the fibred chart 
y\ 2* |J| = ^ s i s < fc) in the domain W k = tt^ (W) C J k (ir). Coordinates z\ in the 
fibers of jet bundles are defined by the condition Zj(j k s) = d 1 s fl (x). 

For k = 1, . . . , oo the space J fc (vr) — > X of k-jet bundle is endowed with the contact (Cartan) 
distribution Ca k defined by the basic contact forms 

uf = dy l - z^dx 1 *, . . . , uj\ = dz} - z} +1 ^dx u , \I\ < k. 

These forms generate the contact ideal C k C Q*(J k (7r)) in the algebra of all exterior forms. A p- 
form is called Z-contact if it belongs to the l-th degree of this ideal (C k ) 1 C Q,*(J k (Tr)). 0-contact 
exterior forms are also called horizontal (or ir^ -horizontal) forms (or, sometimes, semi-basic 
forms). 

Any exterior /-form uj 1 on the space i7°°(7r) has the decomposition oj l = Y2s+r=l u ^ S into its 
s-horizontal, r-vertical components cj( s ' t ' ((s, r)-forms). Denote by ^ s ' r \j°°(Tr)) the space of 
all (s,r)-forms on J°°(7r). 

Bellow we will use the presentation of the exterior differential d as the sum of horizontal 
and vertical differentials dy. for a g-form v G J7*(J fc (7r)) its differential dv lifted into the 
J fc+1 (7r) is presented as the sum of horizontal and contact (vertical) terms: dv = djjv + dyis, 
see PUCE]- Operators da, dy are naturally defined in the space f2*(J°°)(7r)). 

We recall that these operators have the following homology properties 

d\ = d v = dydu + dudy = 0. 

In particular, for a function / G C°° (J°° (it)) (depending on the jet variables z\ up to some 
degree, say |/| < k), 

df = (d^f)dx^ + Y, f,z}^ (2-2) 
|J|>o 

where 

d,J = d x ,f+Y^i + iA}f ( 2 - 3 ) 

|-T|>0 

is the total derivative of the function f(x,y,z) by x^. The series in the formulas (|2.2p . (12 .3D 
contains finite number of terms: 1/1 < k. 
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3 Balance systems, form Kc and the Helmholtz coindition 

Let 7r : Y — > X n be a configurational bundle and G a pseudo-Riemannian metric in X. Balance 
systems of order k is the system of equations for the fields y l (x) of the form 

dpF? + if Ag !(U = IIj, i = 1, . . . ,m. (3.1) 

Here x^, fj, = 0, . . . , n are local coordinates in X - part of the fibred chart (W, cc' 2 , y l ) in V, Ac = 
ln(\/|G|). Densities Ff, fluxes Ff~, A = 1, . . . , n and source/production terms LTj are smooth 
functions on J fc (vr), defined by the constitutive relations C, which determines the properties of 
the physical system (solid, fluid, gas, etc.). 

In the works [131 E] we've associated with a balance systems (|3.ip the n + (n + l)-form 
(modified Cartan form) 

Q c = F^ 1 A + n.a; 4 A 17 + if A n. 

The Euler-Lagrange equations, corresponding to this form Qc 

(j k+1 s)\d@ C = 0, V^G^(J fc (7r)), 

coincide with the balance system (|3.ip . Here d is the "con" -differential acting on a k + (k + 1)- 
form a fc + as follows: d{a k + /3 fc+1 ) = (-da + /3) + d0 (see [13 E]). 

Now we restrict our consideration to the case of a balance system (|3.1|) of order one (i.e. F{~, 
Hi S C 00 (J 1 (7r))) and associate with such system the following (1-contact) (n, l)-form on J 2 (tt) 

K c = (Ftu; l fl + U i u J i )A 11 . (3.2) 

Next we obtain the following version of (local) Helmholtz condition (comp. [T]) for a balance 
system of order one to be Euler-Lagrange system for some Lagrangian L € C°°(J 1 (7r)). 

Proposition 1. For a constitutive relation C of order 1 the following properties are equivalent: 

1) form Kc is closed: dKc = 0, 

2) C is locally Lagrangian constitutive relation: F- 1 = Hi = ^ for the locally defined 
function L £ C°° (J 1 (tt)) . 

Proof. Let the form Kc be closed. Calculate the differential of this form: 

dK c = [(d yj Ik)co j A J + (d yl Ff - d z] Yii)uo l Awj+ (9^)4 A wjj A rj 

= [(dyiUjdyi A df + (d yi F? - d 4 Ui)dy l A tfej + (djF»)dzi A d^] A 77. 

The second line in last formula shows that if we consider Fj~ and Hj as functions of vertical 
variables y 3 , (i.e. if we fix x) then the 1-form of these variables F^dz 1 ^ + Hidy 1 is closed. 
Then, locally, by Poincare lemma, F- 1 dz 1 ^ + Hidy 1 = dyL for some function L(x,y,z). dy here 
is the differential by the variables y J , z 3 v . It is easy to see now that Kc = d(Lrj). 

Proof that the statement 1) follows from 2) is trivial. ■ 

Remark 1. Since dKc = dyKc, the statement of Proposition [TJ equivalent to the local equality 
Kc = dy(Ln), follows from the qualified exactness of the column of the augmented variational 
bicomplex (see bellow), starting at S7 n '°, at the term Q n,1 . 
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4 Balance systems and the variational bicomplex 

Proposition [T] for a balance system (|3.ip to be Lagrangian points at the possibility to interpret 
and study the balance systems using the variational bicomplex, see [TJ [12] . In this section we 
realize this possibility. In presenting the variational bicomplex we will follow pQ. 

The sheaf of fc-contact s + /c-forms on the infinite jet bundle J°°(7r) is denoted by Q, s,k (J°°(7r)) 
or, for shortness, Q s,k . 

We will be using the augmented variational bicomplex - the diagram where all the columns 
and rows are complexes of sheaves of differential forms (right column is the complex of factor- 
sheaves F k = lm(I : fi n ' fc (J°°(vr)) ->■ fi n > s (J°°(7r))). 

o > n ' 3 ■ ■ ■ n n > s — - — > t 3 ► o 



V 



o ' 2 -*s-> n 1 - 2 -2s_> • • • n"- 1 - 2 -2*_> rr- 2 .f 2 ► a^, 

4.1 

d v d v dy d v S v \ 



o ► e ► a - n 1 - . . . n™- 1 - rr° 

Recall also the Euler-Lagrange mapping £: 

5 : fi 71 ' -^ 1 , 5 ([Lt?]) = I(d v Lr)), 

where 7 : f2 ?1 '' s ( J°°(7r)) -4 Q n ' s ( J°°(ir)) is the "interior Euler (homotopy) operator" (see (|4.2p 
below). Mapping 5 is well defined because of the relation dy ° djj = —dH ° dy between vertical 
and horizontal differentials. In terms of this mapping the Euler-Lagrange equations for a section 
s : X — > Y of the configurational bundle n takes the intrinsic form 

(j°°s)itE(Ln) = o, vee^(J 00 (7r)). 

Last row and the last column in the bicomplex (|4.ip form the Euler-Lagrange sequence where 
£ is the Euler-Lagrange mapping 

_> £ _> ^o,o ^ n i,o ^ . . . n n-i,o fc> o n0 ^J^J 2 ^/ 3 ^.... 
We recall (see [HD2]) that the interior Euler {homotopy) operator 

1 : fi n ' s (J°°(vr)) C tt n > s (J°° (vr)) 

that has, as its image, the subbundle J 75 of functional forms of vertical degree s is defined by 
the formula 



I(uS) = -uM 
s 



(j,A),|A|>0 



1 4 

s 



idiW-dpiid , w)+d w/i2 (*a < w) .(4.2) 



Proposition 2. Let iiQ; = H;a/ Ajj + i 7 ^ u 1 ^ An be the K-form corresponding to a constitutive 
relation C of order k, see (|3.2p . Then, the balance system (|3.ip is equivalent to the requirement 
that 

j l s*i^I(K c ) = /or a// £ € ^(J fe (vr)). (4.3) 
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Proof. 



I(K C ) = J A [n iV - d^Ffri)] = & - d^F? - FfXG^u* A V . 

Let £ G <Y(J 1 (7r)). Apply to the last equality and take the pullback by the mapping j°°s, 
s : U —¥ Y being a section of the bundle tt over an open subset U C X. Equation (|4.3|) takes 
the form 

o jh) [d^F? + if A G , M - ILJ o j°° 8 r, = 0. 

Fulfillment of this equality for all (or for many enough) variational vector fields £ is equivalent 
to the statement that s :[/—?> Y is the solution of the balance system (|3.ip . ■ 

Remark 2. The result of Proposition [TJ shows that a balance system (|3.ip is i/ie mosi general 
form of PDEs system that can be obtained from a form of the type (j3.2|) using the infinitesimal 
variational equation (|4.3p in the formalism of variational bicomplex. Similar result for the 
systems of PDEs starting from an arbitrary Lepage form is obtained in |15j . 



4.1 Higher order balance systems 



Balance systems (|3.ip have found their natural place in the variational bicomplex (|4.1|) . Applying 
the mapping / to the form Kc, contracting it with a variation vector field and taking the pullback 
by a 2-jet of a section we restore the balance system (13, ip . 

Thus, it seems natural to inquire, what higher order construction corresponds to the form K 
and which system of PDEs corresponds to this construction. It may lead to the definition of 
a canonical form of higher order system of balance equations. 

Consider an arbitrary (n, l)-form on the J°°(7r) 



K= [ £ifu; A 

|A|>0 



At/, 



(4.4) 



where F t A G C°° (J°° (ir)) . 

Calculate the functional form I(K) assuming that when a = 0, d a is the identity operator. 
We get 



I{K) = lo 1 A 



w 1 A 



£(-1)1^^ K) 



|E|>0 



oj A 



E(-1) |S '^(^) 

|E|>0 



E(-i) |e| ^(^igi)^ 

|E|>0 



Forming the interior derivative of obtained (n, l)-form with a vector field £ € X(J°°(ir)) and 
taking pullback of obtained horizontal n-form by the oo-jet j°°s of a section s : U — > Y we get 
the system of equations for section s 



if *)* E (-^(Od^FplGD^-r, = 0, 
|S|>0 ' ' 



i = 1, ... ,m. 



Role of the source term in these m equations is played by the functions Ff . 
Thus, we have proved the following 
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Proposition 3. For a (n,l)-form K of the form (|4.4p on the infinite jet space J°°(ir) and 
a section s : U —> Y where U C X is an open subset of X, the following statements are 
equivalent 

1) for all (or, what is equivalent, for m linearly independent) vector fields on J°°(7r) such 
that det(y (£,•)) / 0, 

j l s*i^I(K) = 0, 

2) section s is the solution of the system of "higher order balance equations" 
E |G|) =*i, i = l,...,m. (4.5) 



|S|>0 1 



Example 1. Consider the form with F- A = M- for a function L € C°°( J°°(7r)). Then, the 



dz 



A 



dy 

|S|>0 



system of PDEs (j-4. 5f) takes the form 
dL ^ 1 , « / <9L 

4 

Here (— d)s = ( — l) S cfe. For |G| = 1 this system coincide with the higher order Euler-Lagrange 
equations, [3j equation (3.27)] or [12j Section 4.1]. 

5 Vertical splitting of a balance system 
5.1 Vertical homotopy operators 

Recall the construction of the vertical homotopy operator hy S : Q r,s — > f2 r,s_1 defined by I. An- 
derson p]: 



f 1 1 

h V S M = J -K^R^dt, (5.1) 

where R = y l d y i is the Liouville (vertical radial) vector field and 

pr R = tfdyr + + • • • + z\d z i A + ■■■ 

is its prolongation to J°°(7r). The flow of this vector field on J°°(it) is the 1-parameter family 
of diffeomorphisms of J°°(7r) 

$ e (x,y,z) = (z;e e y,e e z). 

It is proved in [T] that for any uj £ Q, r ' s , 

u = d v [h\f(uj)} + h\} s+1 (d v uj). (5.2) 

In the proof of this formula it was shown [TJ equation (4.10)] that the integrand in the for- 
mula (|5.1|) is equal to 

^L^CVR^)! [x,y,z] = t s ~ 2 (i piR uj)[x;ty,tz} = t s ~ l i pr R uj[x;ty,tz}. 

Here u[x;ty,tz] is the form obtained by evaluating the coefficients of the form uj at the point 
[x;ty,tz]. Using this in (|5.ip we get 

hfy s {uj) = \ t s ~ x i wR uj\x;ty,tz\dt. (5.3) 
J o 
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5.2 Decomposition of the form Kc, quasi-Lagrangian 



We will use the formula (|5.3p in the cases r = n, s = 2, 1. 
Let u = K c = IW A r) + F/*w* A?jG ft"- 1 . Then, 

7 prfl a;[x;iy,iz] = IIj(a;,ty,te)ip r ii(a; J A 77) + F?{x,ty,tz)i px A 77) 
= [lL(x, iy, tz)^ + if (x, iy, 77, 

and hy 1 Kc[x;y, z] = Lrj, where the coefficient 

L(x, y, z) = / [y^x, iy, tz)dt + z'if (z, iy, te)] dt (5.4) 

JO 

of the volume form 77 will be called the "quasi-Lagrangian" of the form Kq in the sense of 
decomposition (|5.2p . 

Applying the vertical differential to hy 1 Kc[x;y, z] we find the "Lagrangian part" of the 
form Kc 

d v h n /K c [x- y, z] = jf * (n fe [] + ty% iyk [} + tz^ yk []) dtio k A r, 
+ [ (tym^W + F/[] + A 77 

+ ]T /'Vvn. 4 D+^F.^D)^4Ar ? . (5.5) 

A||A|>1 1/0 V ' A J 

To shorten the formulas we will use the sign [] after function saying that the arguments of 
this function are (x, ty l ,tz 1 ^, . . . , tz^). 

Let 00 = Af^u\ Awj Ai/G O™' 2 be any (n, 2)-form. Then the integrand in (|5.ip is equal to 

ti piR uj[x,ty,tz] = ti pr R [Aff (x,ty,tz)u}\ Aw J E A 77] = L4-- S (x, iy, - 4 w a) A 7?. 

Using this in (|5,3p we get 

h n /uj = (J tAtf[x,ty,tz]dt) (44 - 44) A 77. 

Apply this for u) = dyKc- We have (using the equality dyuj l a = 0), 

d v K c = (d v Ih) A u* A 7/ + (dyFf ) A w* A 7/ = II j 4 A w * A »7 + F M 3 - t^AwlA 77. 

In this formula |S| > 0. Combining this result with the previous calculation we get the comple- 
mental, "pure non-Lagrangian" , term in the decomposition (|5.2p of the form Kc'- 

h n /(d v K c )) = (J m i ^(x,ty,tz)d?j [z^J -yV s ) A 77 

+ (J ^A^tyMdt] (44 - 44) A v 
= U t[( y m hyJ + 41L 4 o + 4n. 4 []) - iyii,,, - ^/-;;; ;/ v (/ /| w * A v 
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+ (J t[-V%, z * A - 4*J, z i)dt\ uj\ A n. (5.6) 

In the last expression [A| > 1 . 
Thus, we get 

Proposition 4. Vertical splitting (|5.2p of a (n, l)-form Kq = Ho/ A n + uj^An has the form 

Kc = Kc,Lag + -Kc.nLag, (5.7) 

where the Lagrangian part -Kc.Lag is given by (|5.5p . lo/izZe i/ie non-Lagrangian part i£c,nLag *s 
gu>en in (15. 6p . 

In the appendix we present several examples of the Lagrangian components (quasi-Lagrangian 
and the corresponding Euler-Lagrange equation) of several well known nonlinear equations 
having the form of balance equations (see Table [T]) . 

5.3 Properties of vertical decomposition 

Consider the bicomplex (|4.ip over a domain W C Y which is vertically star-shaped and covered 
by one fibred chart, see [12] . An important special case of this situation is the case where 
7r : Y — > X is the vector bundle. Then, the vertical complex with the vertical differential dy is 
dy-ex&ct ([121 Theorem 5.58] or [H Chapter 4]): any dy-closed (n, r)-form over W is exact. In 
addition to this, vertical homotopy operator is defined in the domain W°° = ^^q(W). 

Denote by Z n,r the subspace of dy-closed forms in {l n,r . Introduce the subspace Q n,r = 
h n,r+i z n,r+i c fi n,r_ In the decomposition u = h n ' r+1 d v uj + dyh n > r UJ, 

h n ' r+1 d v uj G Q n ' r , dyh n ' r uj e Z n ' r . 

Proposition 5. Locally, in the domain W°° over a vertically star-shaped and covered by one 
fibred chart domain W C Y the following is true: 

1. The decomposition (K -decomposition) 

n n ' r (w°°) = Q n ' r (w°°) e z n ' r (w°°) (k) 

is direct over the domain W. 

2. The presentation of a (n,r)-form uj 

uj = h n > r+1 d v uj + dyh n ' r uj (5.8) 

as the sum of terms in the decomposition is unique. 

3. The mapping V = h n ' r+1 o dy : Q n ' r — > J7 ra,r is the projector: V 2 = V. 

4. The mappings 

Z n ' r+1 Q n ' r and Q n ' r ^ Z n ' r+1 

are linear isomorphisms inverse to one another. 

5. For v G fl n > 1 , h n ' r v G Z n ' r ~ x iff v G Q n ' r . 

6. An (n, l)-form uj is anti- Lagrangian, to = h n,r+1 dyto <3> h n,r uj G TT*Cl n (X). 
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7. Over the domain W°° the following direct decomposition (F -decomposition) at the level of 
functional is valid 

F 1 (W co ) =l(Q n ' 1 ) ®£(n n '°), (F) 

where £(Q n '°) is the space of Lagrangian functional forms. 

Proof. To prove the first statement, it is sufficient to show that the intersection of the two 
subspaces in the decomposition (fK|) is zero. Let oj G Q n,r f] Z n,r (W°°), i.e. 

u = h n ' r+l a = d v X, 

where a = dya 1 ', a' G Q n ' r and A G Q n,r ~ 1 . Applying dy to this decomposition of oj we get 
dyuj = = dyh n ' rJrl a and, therefore, dyh n,r+1 dyo' = 0. Using the decomposition (|5.2|) for a' 
we have h n,r+1 dya' = a' — dyh n,r a' . Substituting this equality into the previous one, we get 

d v h n ' r+1 dya' = dy(a' - dyh n ' r a') = dyo' = 0. 

It follows from this that a' G Z n,r and that a = dya' = 0. Therefore, oj = 0. 

The second statement is another form of the first one. The third statement follows from the 
second one if we apply V to the decomposition (|5.8p . 

To prove 5th statement let v € Q n,r , then, v = h n,r+l dyv + dyh n,r v. By unicity of such 
presentation, dyh n ' r u = 0. This proves the 5th statement and part of statement 4. 

Let a € Z n ' r+1 , then a = d v u, v G Q n > r . We have h n ' r+1 a = h n > r+1 d v v and 

dyh n > r+1 0- = dyh n ' r+l d V V = dy(l- dyh n ' r+1 )v = dyV = a. 

Combining this with the proven statement v = h n,r+l dyv we finish the proof of forth statement. 

Next statement follows from the fact that Ker (dy) in O n, ° coincides with ir*Q, n (X). 

To prove the last statement that the sum (jFj) is direct we use the fact that the homotopy 
operator hy and the horizontal differential dfj anticommute: If u) G Q and if I(oS) = I(dyLrj) for 
some function L, then uj — dyLn = d^fi for some j3 G Q n ,1 . Applying decomposition (jiT)) to the 
form dff/3' daft = h n ' 2 'dy djj ft + dyh n ' 1 dn ft substituting it into the expression u = dyLn + dnP 
for u and using unicity in the decomposition §K§ we see that oj = h n,2 dydH(3 = dnh n ~ 1,2 dy/3. 
As a result, I(oj) = 0. ■ 

Remark 3. I am grateful to the the referee for the remark that if the configurational bundle 
7r : Y — > X is vector bundle, decomposition (|5.2p and (|5.8p are globally defined and depend only 
on the structure of the vector bundle ir. 

Remark 4. Projectors I and V do not commute. It follows from the fact that the forms in the 
image of / are of the type PiOJ 1 A r\ while the forms in the image of V have also terms A^oj\ A n 
with |A > 0. 

Remark 5. It would be interesting to know, if the restriction of hy T to Q n ' r is the monomor- 
phism? Is it onto the Z n ' r ~ 1 ? 

5.4 (n, l)-forms K with zero Lagrangian part 

Here we study the (n, l)-forms Kc which have zero Lagrangian part. 
We will be using the following simple lemma. 
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Lemma 1. Let f(u l ) be a function of m variables u % , analytical in a connected domain D C 
containing the origin 0. Let s G Z be an integer number. If 



\ t s f(tu)dt = c — const (respectively zero) 
Jo 



in an open subset O C D containing the origin, then f(u) = (s + l)c — const (respectively zero) 
for all u £ D. 

Proof. Present function / by its Taylor series centered at the origin 0; 

f(u)= c ^ 7 - ( 5 - 9 ) 

/,|/|>o 

where u = u 1 11 ■ ■ ■ u m %m is the monom corresponding to the ordered multiindex /. This Taylor 
series absolutely converges in some disc ||tt|| < R, where R > 0. Rearrange this series collecting 
terms cju 1 with III = k: 



oo 



/(«)=x;« w , 

k=0 

where vP^ is the homogeneous part of series (|5.9p of degree k. This series also absolutely 
converges in the disc ||u|| < R. We have, now 



{k},k 



i(tu)=j2' 

k=0 

As a result, 

rl oo fl oo 

JIJ k=0 J[> k=0 I,\I\>0 

Last series has the same radius of convergence as the series (15, 9p . Under the condition of 
lemma, series (|5.10p equal constant c (respectively zero) for all u in the disc ||u|| < R. This 
series presents an analytical function and is equal to the constant (respectively zero) in the disc 
< R. Therefore, series (|5. lOj) have all but the zeroth (respectively all) Taylor coefficients 
equal zero. Then, the same is true for the series (|5.9|) for f(u). Being equal constant (being 
identically zero) in an open subset of its domain, function / is the constant (identically zero) in 
the whole connected component of its domain containing zero. This proves the lemma. ■ 

Applying this lemma to the quasi-Lagrangian L given by (|5.4p presented as the sum of 
homogeneous components 

„1 oo 

L= / y2[y i n ik (x,y,z) + zi 1 F< M k (x,y,z)]t k dt 
Jo k=o 

we prove the first statement in the following 

Proposition 6. Let Kc = rij(rE, y l , z^uj 1 Arj + F^(x, y l , %*)&* A 77 be a (n, \)-form analytical by 
vertical variables y l , z 1 ^ in a connected domain D C W n+mn containing the origin. Then, 

1. The quasi-Lagrangian L corresponding to the form Kc is a function of x only (identically 
equals zero) if and only if 

y l Hi(x,y,z) + z^Fj 1 (x , y , z) = cfi(x) (respectively, = 0). (5-11) 

2. If this equality holds, then L = 0. 

Proof. Second statement follows from the first and the fact that the analytical function in the 
left side of (|5.11j) vanish when y % = z 1 ^ = 0. ■ 
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6 Balance systems of zero order 

6.1 Balance systems whose quasi-Lagrangian L has zero Euler-Lagrange form 

Here we look at the case where Lagrangian part does not enter the functional form of Kq (pure 
non-Lagrangian case). 

We consider here only balance systems of order zero (called the RET systems in continuum 
thermodynamics [H]). Thus, the density- fluxes and the source terms IL depend on the 
fields y l and, possibly, on x^. 

Lagrangian functions L whose (n, l)-form dyLAr] is of order zero are linear by the derivatives: 
L = Lq(x, y) + zJ 1 L f ^(x, y) and the corresponding Euler-Lagrange systems are quasilinear systems 
of the first order 

~ L tjv% + ( W " L t^) = °> i = l,---,m. (6.1) 

Consider now an arbitrary form K = Hiuj 1 Arj + Ff^oo 1 A rj of order zero. Quasi-Lagrangian L 
corresponding to this form is 



L= f [y k U k (x,ty) + z k Fj^(x,ty)]dt. 
Jo 



This function depends linearly on the jet variables. 

Applying the Euler operator / to the form dv{Lrf) we get 



dzj, 



In the zero order case, 



= jf 1 {HQ + ty k U Kyl [] + tz k u F^[]} dt - jf 1 {F?J + tz k F» yk []} dt. 

In the homogeneous case, where there is no no explicit dependence of F? on x^, the first term 
in the second integral vanishes and we get 

L„ i -d II .S^r = 



C {HQ + ty k U k ^[]}dt + £ {tz k u F^ yl [] - tz k F» yk []} 
jT {HQ + ty k U k>yi [}} dt + z k J Qf * Fjfldt) t - (jf* Ft[]dt 



The variables z k are independent and are present only as factors in the second term of last line, 
so, L yt — <^ Jpr = if and only if both terms - one containing Ilj and the others containing F? 
all vanish. In particular, we have to have 



ffZUdt) =([ 1 Ft\\dt) , 

JO J , v i \J0 J , v k 



V i, k. 

This condition has the form of mixed derivative test and as a result, Euler-Lagrange equations 
for Lagrangian L vanishes if and only if the following conditions holds: 



J [n^+ty k n ktyi \\)dt = — i J y k U k []dt = 0, 
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locally, for some functions G^(y l ), 



•l 



dC 1 



J F k( t v) dt= -^yki p = l,...,n, k = l,...,m. 

Using Lemma Q] we show that in the analytical case the first condition is equivalent to the 
condition y k Hk(x,y) = c(x) and, due to the analyticity by vertical variables, to the equality 
y k Uk(x,y) = 0. More than this, in the analytical case we can rearrange terms of Taylor series 
for the functions F£ and G ll (y) collecting together monoms y = y lH ■ ■ ■ y m %m with the same 
value of \I\ = k: 

oo 

V 1=0 

Substituting this into the second condition and using Lemma [H we write it in the form 

i=o y p J{) v 

Equating homogeneous terms of the same order in both parts we get 

1 dGfLi 

iTTtw = -sjf 1 ' p = o,i„... 

Multiplying by (p + 1) and taking summation by p we get the formal equality 

OO OO q / oo \ 

K = £ *W = £((?+ v<%*.i)* = A E(p + • 

p=0 p=0 y \p=0 J 

It is easy to see that the series on the right side has the same radius of convergence as the 
corresponding series in (|6.2j) . Denote obtained analytical function by 

oo 

G»(y) = ^2(p + l)GZ +1 (y). 

p=0 

As a result, we have proved the following 

Proposition 7. Let the balance system (|3.ip be analytical and of order zero (RET case). Then 
the Lagrangian part of this balance system in the functional form splitting ([7]) vanishes, or, in 
other terms, I(Kq) £ I(Q n ' 1 ) if and only if the following conditions are fulfilled: 

1) y k Il k (x, y) = c- const, 

2) Godunov condition (see below Definition^: the form F^uj % A rja is vertically closed, i.e. 
locally, F^oj 1 A rjn = dviG^n^) for some analytical functions G^(y l ), fj, = 1,. . . ,n. As 
a result, 

dG^ 

F k^ = ~dyT' /x = l,...,n, k = l,...,m. 



Remark 6. Notice that the Godunov condition has the form dyF^oj 1 A = 0, is covariant 
and, therefore, globally well defined. The second condition: y % Hi = can be written, using the 
Liouville vector field R = y t d y i, in the form i^Kc = 0. Thus, this condition is globally well 
defined on the vector bundles ir : Y — > X where the transition transformations are linear on the 
fibers of n. 
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6.2 Systems of Godunov and the decomposition of the (n, l)-forms 
of order zero 

In 1961, S. Godunov [5] defined an interesting class of balance systems of the form 




where G^{y) are some functions of the dynamical fields y l . These systems are generalizations of 
symmetrical hyperbolic systems of Friedrichs. Later on, S. Godunov and his collaborators used 
these systems as a starting point for constructing Galilean invariant and thermodynamically 
compatible models of a large class of physical systems [H [7] . This was done by adding terms 
depending on the derivatives 2' to the fluxes F^, A = 1,2,3 and introducing the source 
terms Ilj. At the same time, G. Boillat and T. Ruggeri proved that the Godunov systems (j6.3j) 
present the canonical form of a RET (zero order) balance system with the entropy balance if the 
entropy density is convex (condition guaranteeing the thermodynamical stability) |11| Chapter 3, 
Section 2]. 

Thus, it seems natural to introduce 

Definition 1. A balance system (13. ip of order zero and the corresponding (n, l)-form K will 
be called Godunov type system if the density-flux part of the system has the form 

„u 9G" 
dt 

for a functions & e C°°( J fc (vr)). 

Thus, the result of Proposition [7J can be formulated as follows: 

Corollary 1. An analytical balance system of zero order is pure non-Lagrangian if and only if 
it is of Godunov type and its source part Yliu' 1 A ij satisfies to the condition 

y k U k (x,y) = 0. 

Next we project this result to the complex of functional forms (most right column in (|4.ip ). 
where 5y = Iody is the induced vertical differential. Using the vertical splitting (TF|) of functional 
forms (see also [TJ Chapter 4A]) we get the following 

Corollary 2. Let the domain W <ZY satisfies to the conditions of Proposition [SJ Then, any 
analytical functional form K 6 F l of order zero in the domain W is the sum 

I{K) = G k +£ Lk 
of the functional forms generated by: form of Godunov type 

G K = [(IT - L K ^J + (if - ^, 4 K] A n 
and the Euler form 

£l = dyLx A rj 

for the quasi- Lagrangian Lk- This decomposition is unique {in the domain W). 

Proof. Present the (n, l)-form K as the sum K = K\ + dv(Lxr]) as in (|5.7p . Component 
K\ € Q n ' 1 in this representation has the property that Lk 1 € C°°{X) and, therefore, £jh k = 0. 
From the previous corollary it follows that K\ is of Godunov type. Then, since L 2 = L, K = 
L(K) = L(K\) + £^ . Now we rename L{K\) = Gk and the proof is complete. ■ 
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Comparing systems of equations corresponding to the Lagrangian component dv{Lrj) and 
the Godunov component of a (n, l)-form K = Tlico 1 A rj + F^oj 1 ^ A r/ with IT, € C°°(Y), we 
see that the tensor of the principal parts A^(x,y)dyj (of i-th equation) is symmetrical 

(by ij) for the Godunov component of the balance system and anti-symmetrical (see (|6.ip ) for 
the Lagrangian component. 

Let now K = Hi(x,y)u} 1 A i] + F/^x, y)u;^ A r/ be a (re, l)-form of order zero of the type 
considered here. Let L = Lo(x,y) + ztL^{x,y) be the corresponding quasi-Lagrangian. Then 
L\ = f Q Ff(x,ty)dt. Godunov component Gk of the form K has the functions Ff = Ff — 
Ak>* = Ft ~ Jo F^(x, ty)dt as its densities (/j, = 0) and fluxes (/j, = 1, 2, 3). Matrices 

Ff , = F. M , - / Ax,ty)dt 

are symmetric by (see Definition [1]) . As a result, the system of balance equations corre- 
sponding to the form Gk is the first order system of PDEs M^(x, y)y 3 ^ = bi(x, y) with symmet- 
rical matrices M^- . If the matrix Mf- = F? j is positive definite, this system is the symmetrical 
hyperbolic by Friedrichs [11]. Thus, we get 

Corollary 3. Let K = Ilj(x,y)a; 1 A rj + Fi(x,y)uj l „ A r/ be an analytical by y l (re, l)-form of 
order zero and let the symmetrical matrix = F® yj is positive definite. Then the system of 
balance equations corresponding to the Godunov component Gk of the form K is the symmetrical 
hyperbolic system by Friedrichs. 

We finish this section with two examples of calculation of quasi-Lagrangian for two balance 
systems - the system of elasticity equations for a hyperalastic body and the reduced system of 
equations of 2-dimensional ideal plasticity. 

6.3 Example: hyperelasticity 

Consider the system of equations of hyperelasticity JTH} Section 1.6]. 

Let M 3 be a material manifold - a connected open 3-dimensional manifold endowed with 
the (reference) Riemannian metric go- Local coordinates in M will be denoted by X , / = 
1, 2, 3, time t by X°. A deformation history of the body M is a time dependent diffeomorphic 
embedding <f> : M — > E 3 of M into the Euclidian 3-dimensional (physical) space E 3 endowed 
with the metric h. Local coordinates in E 3 will be denoted by x l , i = 1,2,3. As the dynamical 
variables of Elasticity we choose the velocity vector field over (f>: v = and the deformation 

gradient F, F\ = J^- - (1, l)-tensor over the mapping (p. 

Constitutive properties of the body are determined by the strain energy function id = w{F 1 ^). 
Strain energy w depends on the material metric go and the physical Euclidian metric h |10| . 

Elasticity equations have the form of the balance system of order zero, where the first equation 
is the linear momentum balance law and the second one - compatibility condition, 

d x oVi + d XI (^j = bi(X), d x0 F} - d XI {5lv k ) = 0. (6.4) 

Dynamical fields are y l = v l , i = 1,2,3 and y aA = F%, a,A = 1,2,3. Jet variables will be 
denoted by zi, i = 1, 2, 3, fi = 0, 1, 2, 3, corresponding to the derivatives of velocity components 
and z^ A , a, A = 1,2,3, fi = 0,1,2,3, corresponding to the derivatives of components of the 
deformation gradient. Basic contact forms in the second jet bundle J 2 (ft) are 

J = dv l - zldX u , u aA = dF% - z a v A dX v , 
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i = dv^ - zj^dX" , u aA = dF%,» ~ ZA^dX" 



The form Kc in this case is 

K c = (W A t?) + (vn4 + w, F1 LU aA ) A ?? + {F%ujI A - 5%v A u; B A ) A rj. 

To calculate the quasi-Lagrangian L = J [y l Ui\\ + ZqF^\\ + z^-F/^dt we write 

[\y i U i ]dt = v i b i , 
Jo 

which has the meaning of mechanical power of the bulk force bidx 1 (we assume that bj(X) are 
functions of X only). 

Next, using the compatibility equation, we get for the terms corresponding to the elasticity 
equations 

J (4*?U + AF^dt = J X (4t Vi + z\w^ A [])dt ~ Ividtv* + v\ xA J w^Jdt 
= O x o O^f-) + (d x oF\)d Fk J r'wWdt = d x o (^f- + jf * t^wWdt) . 

We recognize the time derivative of the energy density of elastic body (sum of kinetic energy 
and strain energy) . Finally, using F B = v a x B , we get 

f Q [4 B F°aB + z a A B F* B ]dt = jf* [zfF%t + z a A B S A vH]dt ~ \[F%pF% + F a B xB v a ) 

= \{v: B F a B + v a Fl xB ) = l -d xB {v a F a B ). 

Combining we get (partial derivatives here coincide with the total derivatives) 

i2 r l 



L = v% + d x o ^M- + t^wftdtj + d xB ^\v a F^j 



Thus, the quasi-Lagrangian is the sum of the trivial Lagrangian and the source term that 
generates the right side in the elasticity equation (|6.4p . Notice that the trivial part of L has the 
form of the left side of a balance law similar to the energy balance for the elastic body. 

6.4 Example: 2-dimensional ideal plasticity 

System of equations described the plane strained state of the medium with von Mises conditions 
can be transformed, by change of variables (see to the system of balance equations 

U£ + -v = 0, v tV + -u = (6.5) 

for the functions y 1 = u, y 2 = v of variables x° = £, x 1 = rj. For this balance system, F® = u, 
Fl = 0, F 2 ° = 0, F\ = v, n x = -\v, n 2 = -\u. 
The quasi-Lagrangian of this system is 



2 



uv. 



Thus, in this case the quasi-Lagrangian is also sum of total divergence and of the expression 
producing the source terms in the balance system (|6. 5[) . 
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7 Conclusion 

This works shows that the general systems of m balance equations (|oM|) for m dynamical fields 
y l represent the natural type of systems of PDEs appearing in the framework of variational 
bicomplex. Any such system is associated with the (n, l)-form Kq = ILu;* A rj + F^cu 1 ^ A r\ € 
Q n ' (J°°(7r)). It is proven that the vertical homotopy formula delivers, in appropriate domains 
or, for vector bundles, globally, the unique representation of Kc as the sum of the Lagrangian 
term and the "pure non-Lagrangian" form. The corresponding splitting of the functional forms 
follows from the decomposition of (n, l)-forms. It is shown that for zero order forms Kq, whose 
quasi-Lagrangian L is trivial (I(dvLrj) = 0), the functional form I(Kc) is of the type introduced 
by S. Godunov in 1961 [5], and, later on, identified as the canonical (symmetrical hyperbolic) 
form of balance system in the Rational Irreducible Thermodynamics, [llj . 

In the continuation of this work we will study the decomposition (|5,7p of the functional forms 
into the Lagrangian and non-Lagrangian parts for the balance systems of order one. Another 
natural direction of future work is to study the question (suggested by the referee) : is it possible 
to obtain a global version of the splitting ([71) ? Finally, it would be interesting to extend results 
obtained for the case of analytical constitutive relations (density /flux and source terms Fj 1 , 11,) 
to the smooth (C°°) constitutive relations. 

A Examples of vertical decomposition 

In this appendix we present a few examples of the Lagrangian component of several well known 
nonlinear equations having the form of balance equations. As these examples show, the terms 
in second order spacial derivatives arrive from the Lagrangian component. The source terms in 
examples of Subsections 16.31 16.41 also come from the Lagrangian component. The same is true 
for the source term in the last equation below. It is possible that observations of this kind will 
be useful for the construction of thermodynamically compatible balance systems for different 
physical systems (comp. [HE]). 



Table 1. 



Name 


Equation 


quasi-Lagrangian L 


Lagr. part of eq. 


conser- 
vation law 


d t u + d x C(u) = 


dk (4) + d * (Jo Jo C{ty)dydt) 


= 


KdV 
equation 


d t u + d x (3u 2 + u xx ) = 






= 


Burgers 
equation 


d t u -d x (^- + u x }j = 




u xx = 


Filtration 
equation [5] 


(u - Au),< - (u >xA ) tX A = 






Ati = 


cubic Schro- 
dinger equation 


iu,t + + a\u \ 2 u = 


d t {iu 2 ) + i||Vu|| 2 -fuM 2 


Au+ f \u\ 2 u = 
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